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The principle of detailed balance is at the basis of equilibrium physics and is equivalent to the 
Kubo-Martin-Schwinger (KMS) condition (under quite general assumptions). In the present paper 
we prove that a large class of open quantum systems satisfies a dynamical generalization of the de- 
tailed balance condition {dynamical detailed balance) expressing the fact that all the micro-currents, 
associated to the Bohr frequencies are constant. The usual (equilibrium) detailed balance condition 
is characterized by the property that this constant is identically zero. From this we deduce a simple 
. . . and experimentally measurable relation expressing the microcurrent associated to a transition be- 

04 , tween two levels em en as a linear combination of the occupation probabilities of the two levels, 

' with coefficients given by the generalized susceptivities (transport coefficients). Finally, using a 

' master equation characterization of the dynamical detailed balance condition, we show that this 

condition is equivalent to a "local" generalization of the usual KMS condition. 

PACS numbers: 05.60.Gg, 05.30.-d,03.65.Yz, 
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I. INTRODUCTION 

> . 

QQ , To understand non-equilibrium phenomena is one of the most important challenges of modern physics. The mono- 
00 ' graphs |^ summarize the early developments in this direction and, after them, several endeavors were made by 
many authors to construct a satisfactory description of non-equilibrium phenomena from the stand point of micro- 
scopic physics (cf. e.g. ||, ^, ^, ^, |l^). As pointed out by many authors (see for example (l^), the most crucial 
difhculty of the problem is that we lack a good characterization of non-equilibrium states whereas we have criteria 
for the equilibrium case: detailed balance, the KMS condition, stability and so on. In the present paper, starting 
from some physically interesting situations we deduce a general characterization of a class of stationary states which 
satisfy a condition (dynamical detailed balance) generalizing the usual detailed balance and KMS conditions. For 
this purpose, we apply the stochastic limit technique |l^ to some concrete and widely studied models 

' and show that this leads to a natural generalization of both the detailed balance and the KMS conditions which 
^ . characterizes a rather wide and interesting class of non-equilibrium stationary states. 

^ ' The first basic idea of the present paper can be described as follows. The most commonly used states in quantum 
Q^, field theory are the Fock (vacuum) or Gibbs (equilibrium) states. When a field in such a state interacts with a 
' discrete system (e.g. an atom) in the stochastic limit one obtains a master equation for the system whose stationary 
• w^ \ state is the ground state of the atom, if the field was originally in the Fock state; while it is the Gibbs state of 
' the system at inverse temperature /3, if the field was originally in its equilibrium state at inverse temperature (3. 
JH ] The systematic development of the theory of stochastic limit (see below and |l^), has revealed that the above 
- - 1 described phenomenon is quite universal namely: for a large class of states (including many concrete examples which 
are neither Fock nor equilibrium) the stochastic limit procedure allows to deduce master equations whose associated 
Markov semigroups drive the system to a stationary state Poo in the sense that, independently of the initial state po, 
one has 



o 

:^ 

Oh: 



hm P'po = Po 



{Pi is the Markov semigroup acting on density matrices). This fact suggests to give a dynamical characterization of 
ground (or equilibrium) states of the system (atom) in terms of their response to an interaction with the environment 
(field) in the stochastic limit regime. This extends to the non-equilibrium regime the approach of In fact, 

from the above considerations it is natural to expect that the analysis of the stationary states of master equations 
associated via stochastic limit to non equilibrium states of the environment, will lead to a new class of states, of 
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discrete quantum systems, which should play for non-equilibrium phenomena, a role analogue to that played by Gibbs 
states for equilibrium phenomena. In the present paper, we prove that this is indeed the case. 

The second basic idea of the present paper is to exploit the main advantage of stochastic limit with respect to the 
old Markovian approximation namely: the field degrees of freedom are not traced away, but they survive in the limit 
as "quantum noise" (or master field). In particular, as shown in [|l^, |l5j, the slow degrees of freedom of the field 
(e.g. the functions of the free energy of the field) survive in the stochastic limit. This allows us to define the energy 
currents in a natural way and to study their dynamics, thus going far beyond the Markov approximation where one 
only obtains the master equation for system observables and looses any control on the limits of field observables. 

We will illustrate our ideas with two models: one is very well studied in the literature and consists of a system 
interacting with two equilibrium thermal reservoirs at different temperatures. The master equation approach to this 
model was discussed in (Q. As already mentioned, this technique cannot be applied to the problem studied in the 
present paper, i.e. the dynamical study of the currents associated to the field because the field degrees of freedom are 



traced away from the beginning. The second class of models is more general (see Sec. VI), because the field, with which 
the system interacts is not in a usual equilibrium state, but in a new class of states in which, roughly speaking, each 
frequency is at local equilibrium at its own (frequency dependent) temperature. Although states of this type have been 
considered in studies of molecular kinetics, we do not know if these states have been experimentally realized. However 
their structure, characterized by local equilibrium at energy dependent temperatures, is a natural modification of the 



usual equilibrium states (see Sec. VII below) and we are confident that the inventiveness of experimentalists is rich 
enough to allow their realization. 

We briefiy describe the general scheme of the stochastic limit technique for Hamiltonians of the form 

H^^^ =Ha + XHi (1) 

where A is a real parameter, Hq is the free Hamiltonian and Hi is the interaction Hamiltonian (see the concrete 



example in the next section). The general idea of the stochastic limit approach |11, |l2| is to introduce the time 
rescaling 

t^t/X'' (2) 

in the solution 

^ ^^tHo^-^tH(^^ (3) 

of the Schrodinger equation in interaction picture associated to the Hamiltonian H^^\ i.e. 

^ U^^^ = -tXHi{t)U^^\ Hi{t) = e''""Hie-''"'>. (4) 



dt 

The rescaling (||) gives the rescaled equation 

Pl/l = -{Hi{t/X^)U^l (5) 

and the limit A — )■ (which is equivalent to A ^ and t oo under the condition that X^t tends to a constant) 
captures the dominating contributions to the dynamics, which, under appropriate assumptions on the model is 
shown to converge to the solution of 

^Ut = -ihtUt, ht = lim \Hi{t/X^) , U{0) = 1 (6) 
dt A^o A 

Similarly one obtains the limit of the Heisenberg evolution 

lim^ Xl^^ := lim U^^J XU^^^ - UjxUt (7) 



where Ut is the solution of (|6|) and X is an observable belonging to a certain class (slow observables, cf. Sec. HI below 
andll).^ 

The main result of this theory is that the time rescaling induces a rescaling 

Qk — > -e ^ 'ok (8) 
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of the quantum field, defining the Hamiltonian (1), which in the present paper will be assumed to be a scalar boson 
field: {[0^,(1^'] = S{k — k')) (the meaning of Lo{k) and 57 will be described in next chapter) and, in the limit A — *■ 0, 
the rescaled field becomes a quantum white noise (or master field) 6n(t, k) satisfying the commutation relations 

[5n(t, fc), bl,{t', k')] = 5nM''^TiS{t - t')5{k - k')5{iu{k) - Q). (9) 

Moreover, if the initial state of the field is a mean zero gauge invariant Gaussian state p/(0) with correlations: 

{alak') = N{k)S{k-k') (10) 

then the state of the limit white noise will be of the same type with correlations 

{blit, k)ba' it', k')) = <5o,n'27r(5(i - t')5{k - k')5{uj{k) - n)N{k) (11a) 

{bnit, k)bl,{t', k')) = 5n,n'2TT&{t - t')5{k - k')5{uj{k) - Q){N{k) + 1). (lib) 

It is now well understood that this scheme plays an important role in the analysis of the limit (^ when X is a 
system operator. In Sec. HI, we describe a new development of the stochastic limit which allows to extend this scheme 
to a class of observables describing the slow degrees of freedom of the field. 

The remaining part of this paper is arranged as follows: In Sec. || we consider a model which drives the system to 
a non-equilibrium stationary state. It describes a quantum system put between two reservoirs at different tempera- 
tures. By analysis of the reduced density matrix with stochastic limit, we show that this system has non-equilibrium 



stationary state which doesn't satisfy the detailed balance condition. In Sec. Ill, we apply the stochastic limit to the 
slow degrees of freedom of the field. This allows to define the currents associated to these degrees of freedom and to 
discuss their properties. In terms of these currents, we define the dynamical detailed balance condition which is a 
generalization of the usual detailed balance condition. In addition, we show that in the linear approximation these 



currents satisfy the Onsager reciprocal relations [ |18| . In Sec. IV, we investigate a master equation characterization of 
this dynamical detailed balance condition, which corresponds to the well-known fact that the usual detailed balance 
condition is characterized by the master equation which drives the state to equilibrium [ p^ , pT| . Then in the next 
section, we introduce the local KMS condition and prove that it is equivalent to the dynamical detailed balance 
condition for the state. In addition, we consider another model in which the system interacts with an environment 
whose state is non-equilibrium and satisfies the local KMS condition. We show that such states of the environment 
drive the system to a non-equilibrium state satisfying the local KMS condition with a non linear temperature function 



which is uniquely determined by the state of the field. Finally in Sec. VII, we summarize the contents of this paper 
and discuss related topics. 

II. DEDUCTION OF THE STOCHASTIC SCHRODINGER, LANGEVIN AND MASTER EQUATION 

In this section, we consider a model in which the system is driven to a non-equilibrium stationary state by its 
interaction with two non-equilibrium boson fields. This interaction is described by the Hamiltonian 

H = Hq + X Hj^, (A is a coupling constant.) (12a) 



H, = / dk 



(g,{k)D,a]^^g*{k)D]a,.k) (12c) 



where Dj and Z3j are operators on the system space, aj^k and a}- j, are the annihilation and creation operators of the 
j-th field (j=l,2) and gj{k) is a form factor. 
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The initial state of each field is a Gibbs state at temperature f3j ^ and chemical potential fij with respect to the free 
Hamiltonian (throughout the present paper we assume LOj{k) — fij > for all k as usual), i.e. the mean zero gauge 
invariant Gaussian state with correlations: 

{0'lk"'j'M)^SjfN{k;Pj,fij)S{k-k'), N{k■pJ,^lJ)^-^-^^J-^^—^—^ (13) 
The Schrodinger equation in the interaction picture is 



where 



Hi{t) = ^''""Hi, 



-itHo 



3 = 1,2 



= E EE / * {9r.irn{k)EUlrn)ay+^^'^^^'^^--^'+gl^^^^^^^^ (15a) 

j=l,2wGF Im •' 

gr.im{k) ^ gj{k){ei\D\e„i), E^{lm) ^ ^ (e^ - cj|ei)(em|er)|er - w)(er|, (15b) 

F = {oj = Ej. - er';er,er' € Spec. (i^s)}, ^I, — {e^' G Spec.(i75); e,.' — € Species} (15c) 

In the following, for simplicity, we assume that Hs is generic, i.e. 

1) the spectrum Space Hg is not degenerate 

2) For any uj \ F^\ = 1, i.e. there exist a unique pair of energy levels e;, e„i € Spec.(i7s) such that a; = — e; 



In such case, ( |l5a| ) becomes 

^^(*)= HY. (.9,,.(fc)£;.at^e+'(--W--)*+5;:.(A;)i?ia,,,e-(--«--)*) (16a) 

where 

93:uj{k) ^ g]{k){ei\Dj\e„i), ^ \ei) {eyn\: for e;, s.t. e„ - e; = w. (16b) 

Giving an Hamiltonian such as ( [l^ ) the stochastic limit technique proceeds in four steps: 

1. Write the associated white noise Hamiltonian (WNH) equation ([T^). 

2. The causally normally ordered form of the WNH equation gives the Stochastic Schrodinger (SS) equation (|is|). 

3. From the SS one deduces the Langevin equation (e.g. ( |2l| ) and (|3l|)). 

4. Partial trace of the Langevin gives the master equation (e.g. (p2|)). 

In the following, we shall describe the results of these steps for our models and we refer to jl^ for a detailed description 
of the steps necessary to achieve these results. 

Applying stochastic limit as explained in Sec. ^ we obtain the white noise Hamiltonian equation 

^Ut i^^^k^ + Elh-.jJ Ut (17a) 



where 



dt 

3 = 1.2 uieF 



bt:,,u = / dk g*.,Jk)bt..j,Uk), bf.j.Uk) = Hm ie-^(".- )*/^'a,- fe. (17b) 
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Notice that the state of the hmit white noise will be of the same type as ( jl^ ) but with correlations 

H.jjk)br..j',^ik')) = 6j,,6^^^27r6it~t')d{k~k')S{io,{k)-Lu)Nik;P,,^,,) (18a) 
{bt:jAf^)bl^y^^,{k')) = 6J,,6^^^2^^6it-t')6ik-k')S{io,ik)-Lu)iN{k;^3J,^i,) + l). (18b) 

The SS equation associated to the WNH equation (|l^) is 



where 



dUt = -iYl {^^d^k^ + EldBt;,,^ - i (7_ J- ^S^^t + ^X^^^^ElE^) dtj Uf (19a) 

3 = 1,2 



/t+dt pt+dt 
br:,,^dT, dBl^^^ = bl^^T (19b) 

are stochastic differentials and satisfy the Ito table 

dBt;j^ujdBlj,^^, — 2Sjj'Sujuj'Rej-^j^ujdt, dB\.j^dBt;j'^uj' = 2Sjj'dujuj>Rej+j^^dt (19c) 

dtdBf.j,^ = dBt..j,^dBt..f,^' = dBl^^Bl^,^, = dtdBl^^^ = 0. (19d) 
The main physical information is contained in the generalized susceptivities (or transport coefficients): 

7_j> = / dk\gj^^{k)f 



.^-iiN{k;Pj,fXj) + l) 



Lu — U!j{k) — iO 

r |g^.^(fc)|2 gP^(^Ak)-H) 



= t: [ dk |5 (fc)|2-__i ^ Siuj,ik)~u;)~iP.P [ dk ^ , ^ . (20b) 

For an operator X of the system space "Hs, from the SS equation (|l^), one obtains the Langevin equation 
d {ujXUt) = E E {^t [E^ : X] UtdBl^^ + U} [El , X] UtdBt.,,^ 

j = l,2 ueF 

-lm-f^^,^^U}[E^ElX]Ut dt + lmj+^j,^u}[ElE^,X]Ut dtj 

-U} (Re7_j, ^ {{E^El,X} - 2E^XEI) + Re7+,,- ^ {{EIE^,X} - 2EIXE^)) Ut dt^ . (21) 

The Langevin equation with the state for some operators of the field degrees of freedom will be discussed in the 
following section (see (|3T|)). Taking the partial expectation value of both sides of this Langevin equation with the 
state (|l3|), the master equation for the reduced density matrix is obtained: 

- E r-,- (\ {ElE^,Ps{t)} - E^psml) - E (\ {Eu^Elpsit)} - Elps{t)E^^ (22a) 



A = ^ E E {^Hi-,3MEu~IHi+,,^)E^EI) (22b) 
Fzp,^ ^ 2Re E 7tj,'^ > 0, (r=p,^ = for uj <0). (22c) 
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The generator of (22a) has the standard GKSL form |T^. For the ofF-diagonal matrix elements Pmn{t) = {^rn\psit)\^n) 
(m ^ n) we obtain 

at 



ijTin — (6'-,e^-e, — 9-^e„-ei — ^+,ei-c„ + £,-£„) , ^T,^ ~ ^^ItJ,'^ (23b) 

I 

G„„ = ^ (r„/ + r„i) > 0, where = | 



^e„,-ei for Cm > Q 

.ei-t^ for em < e; 



(23c) 



A„upu, Ami^\ . (24b) 



which shows that these elements vanish at i — > cxd whenever Gmn 7^ 0, m,n). 

The diagonal matrix elements pmm{t) describe a classical birth and death process characterized by the equation 

^Pmrnit) = -^((r_,e„-e, + r+,,, )p„„ (t) - (r_,,,_,„ + T +^,^^,,) pii{t)) 
I 

= - ^ (r,™; p,„m(<) - Tim pii{t)) (24a) 

[]; r.,„/ for / = TO 
-Tim for I 

L 

TT- = (for em > e;), or — = (for e, > e™). (24c) 

Notice that this quotient is universal in the sense that it does not depend on gj whenever in the interaction ( p^ ) 
form factors gi do not depend on j {gj — g). When the matrix A has a non-trivial eigenvector associated to the 
eigenvalue, a stationary state exists. In addition, the convergence to the stationary state from any initial state ps(0) is 
guaranteed under quite general conditions (cf. |l^). Notice, that the stationary solution of ( |2^ ) satisfies the detailed 
balanced condition, i.e. 

Pmm ^Im /r>rN 

(25) 



if and only if the coefficients Tmi satisfy 



pu ^ml 



^ Im -L km ^ Ik 



In the non-equilibrium case (^6|) is not satisfied. With this model given by ( |12| ) and (|13|), (^6|) is satisfied only in 
some special cases (for example when both fields have the same temperature and chemical potential, or the system 
has only one Bohr frequency). 

In general, the stationary state of the master equation ( p^ can be described by the nonlinear temperature function 

/3s(em) = — log Pmm. > (27) 
^m 

as 

PS = —^ , Z = trs(e-^-(^-)^-). (28) 

The state is Gibbs state the function (5s{Hs) becomes constant. This fact actually leads to the idea that a rather wide 
class of non-equilibrium stationary states can be treated with such generalized temperature functions. This notion is 



valid not only for the system but also for the state of the field. Indeed, in the Sec. VI below, we will consider another 
model in which the system is driven to a non-equilibrium stationary state by an interaction with a non-equilibrium 
field described by a generalized temperature function. 
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III. MICROSCOPIC CURRENTS AND DYNAMICAL DETAILED BALANCE 



In the previous section, we have investigated the dynamics of a system interacting with fields in a non-equihbrium 
situation and we have already remarked some important difference from the equilibrium case. However one can see a 
more direct and crucial difference through the study of the dynamics of the field degrees of freedom. 



A. Slow degrees of freedom and micro-current 



In order to investigate the dynamics of the field, it is important to notice that some operators of the field degrees 
of freedom, i.e. the slow degrees, survive even after stochastic limit. As we explained in the introduction, the rescaled 
field operators Ofc and a[, become white noise operators denoted by b^{t, k) and blj{t, k) whose commutation relation 
is given by (||). Due to this fact we can intuitively say that the fast degrees of the field become noise (singular) in the 
stochastic limit. However we can describe the time evolution of some of operators of the field in terms of the rescaled 
time even after stochastic limit, and this approach gives us meaningful information on the original dynamics as well 
as on the system operator. Since the stochastic limit is an asymptotic theory, mathematically we have to prove the 
convergence of the dynamics and this has been done elsewhere [12j. In the present paper, we apply the theory to the 
number operator in the model and discuss its physical meaning. 

Let us sketch how to compute the time evolution of the number operator Uk = aj^Ofe under the white noise equation 



^Ut = -^Y. (^-^1 + ^-^*) (29) 

We will illustrate the calculation only in the simplest (Fock) case. The more general states ( p^ ) can be reduced 
to a linear combination of two independent Fock representations (cf. section 2.18). The key formula to apply 
stochastic limit to the number operator of the field — aj^Uk is 

[b^{t, k), uk'] = b^t, k)5{k ~ k'). (30) 

The Heisenberg evolution of n^, after the stochastic limit is described by the Langevin equation 

^ (u^UkUt) = t Ul[E^bl^ + Elbt,^, Ukpt 

= E {^t {EMt, k) - ElbUt, k)) Ut) 

= E (^-(^' kWlE^Ut - UjEMUt, k) + [U},blit, k)]E^Ut - U}El[b^{t, k), Ut_ 
J2 (bUt,k)ulE^Ut - UlEMUt.k) 

dEF 

(jUk)U}ElEM + i^{k)ulElE^Ut)) 5{u{k) - to) (31) 
where in the Fock case 

j^{k)=7r\g^{k)f. (32) 

Taking partial trace ( • ) over the initial state of the system and noise we obtain the evolution equation of the mean 
number of quanta 

j^iUlukUt) = 2 l^mulElE^Ut)5{Lo{k) - u;). (33) 
uieF 

This can be expressed in terms of the time evolution (under the master equation (p^)) of the reduced density matrix 
ps{t), i.e. 

j^iUlukUt) = 2Y5{Lo{k)~-u;)^^{k)ivs[ElE^ps{t)) 

ujEF 

J2 5{Luik)-{e„,-en))trs[juik)\ern){em\ps{t)) (34) 



ujeF 

-i 

ujEF 
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In the case of a general initial state described by ([LOD, the computation is similar, and one get 

j^iUlnuUt) = 2Y,5{Lo(k)-Lo)tTs{{n-Ak)ElE^-l+Ak)E^El)ps{t)) 



2 J2 ^(^(k) - ( l)ps(i)) (35a) 



where now instead of 7t^(fc) given by (32) one has: 

7_^^(fc) = 7-(fc)(A^(fc) + 1), l+.Ak) = lUk)N{k). (35b) 

As the consequence of (p^, once we obtain the time evolution of psit) by solving the master equation discussed in 
the previous section, we find the time evolution of the number operator of field degrees of freedom. 

In order to apply (jsj) to the model discussed in the previous section, we can consider the number operator for each 
field. Defining Uj^k = a] k'^j.k, we obtain 

^^{Uln,,kUt) = 2 ^ S{u;,{k) - w)trs {{j^,^.,{k)ElE^ - j+^^,,{k)E^El) ps{t)) 



uieF 



2 ^ 5{ujj{k) - {em - e„))tr5 ((7-,e„-e„ j(fc)|em)(em| - 7+,e„_£„ j |e„) (en|) /3s(i)) 

51 ^i^j(k) - (fm - en))(7-,e,„-c„,i(fc)Pmm(i) " 1+ ,e^-e„ ,j Pnn{t)j (36a) 



em>€ 

2 



where 



7_,^^,(fc) = 7r|g,,^(fc)|2 (Ar(fc;/3^.,^^.) + 1), 7+^^,,(fc) = 7r|.g,-^(fc)|2 iV(fc; /3„ ^i,). (36b) 



This time dependence of the slow of degrees of freedom of the field is due to the interaction with the system and is 
a direct evidence of the existence of a family of currents passing through the system: one for each proper frequency 
= Em — £n > 0. To investigate these currents, let us define, for each e,„ > e„, the region flmn in fc-space, resonating 
with the frequency LUmn ■= £m — £n which includes all kmn such that 

i^ihnn) - (Cm - £«) = 0. (37) 

Then define the microscopic number current, associated to the frequency ujmn by: 



= 2 (Re7_j,<:„-c„ Pmmit) - Re7+j,<:„_<:„ p„«(i)) , 
Jj,mn '■ — 2 (Re7_j^e^_c„ Pmm " ^Gl+,j,em-e„ Pnn) , (in Stationary state of the system) 

and similarly the microscopic energy current 

Jf,mn{t)-- = ^(/^ dk u;,{k){u}n,^^Ut)) 

= 2 (e„i - e„) (Re7_j,e^_<:^ Pm,n{t) - Re7+j, Pnn{t)) 
Jf,mn ■ ^ 2 (e™ - e„) (Re7_ pmm - R.e7+j\e,„-e„ p„„) , (in stationary state of the system) 

= 2(e„ - e„)7,,„.„ P^n. ,^ ^ ^ 1 - g-t^— ^^'^^ ^ (38b) 

where 

7j,mn = 7^ / rffc l5j,e™-£„(fc)l^('^i(fc) - (e™ - £«))• (38c) 
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The term microscopic here refers to the fact that we define one current for each atomic frequency. We see, from (|3 
that in the stationary state for the system 

Ps{t) = PS 

we have a constant flow of quanta from the field to the system. 

The sum, over all m and n, of our micro-currents gives two macro-currents which coincide with those defined by 
H. Spohn and J. L. Lebowitz in terms of the master equation In fact, as seen in (^8|), these currents can be 
represented with the matrix elements of the reduced density matrix and the generators of master equation like they 
defined (cf. also the formulas (41) and ( ^ ) bellow). However the micro-currents are essential to define dynamical 
detailed balance and the fact that we started from the dynamics of the fields and deduced them gives a physical 
interpretation to these currents. 

Moreover our approach shows that in fact a much stronger condition is satisfied namely: for each Bohr frequency 
Lu F the mean micro-current relative to the frequency lu = — £n is constant. This means that, for each uj d F, the 
flow of quanta from the modes of the field resonating with the frequency co (in the sense of condition (|3^)) is constant. 
Thus the current of quanta in the field is split into a family of independent microscopic currents^ one for each Bohr 
frequency w. In the stationary state each of these microscopic currents is constant: we shall call this fact dynamical 
detailed balance. This condition gives a simple and experimentally measurable relation expressing the microcurrent 
associated to a transition between two levels as a linear combination of the occupation probabilities of the 

two levels, with coefficients given by the generalized susceptivities (transport coefficients). 

The usual (equilibrium) detailed balance condition is the particular case of the dynamical one corresponding to the 
case in which all the microscopic currents are zero. In fact in this case equation (|38|) is reduced to 



tr(|n)( 



n\PS 



Pmm tr( |m)(r7i|p5 



Vj = l,2 



for any k„in satisfying condition (37). From this, by standard arguments, it follows that there exists a constant /3 > 
such that 



Pi = P2 = f3, P-n 



-pan 



'Per, 



so that PS is the Gibbs distribution. 

In the general case the dynamical detailed balance condition is 



This gives, for m > 



2 (R,e^_^g^_e^^j" Pmm -f^*^T-t-,eTn — en ,i Pnn) — Jj^mn 
2Re7_,e^_eo J Pmm — 2Re7+,e„_eo,j poO + Jj,ma 

Jj^mO 



(39) 



^ _ Re7+,e„-eoj- 

Pmm — -fZ Poo 



Re7- 



2Re7^ 



Replacing this into (B9h we find 



2Re7_ 



Re7+, 



gm-eo,J 



Poo 



j.mO 



Re7-,<:n.-eo,i 

2Re7+,,„„e„j 



2Re7_,<:^_eoj. 



Poo 



or equivalently 



7—2 



_Re7_,e„_eoj 2Re7_,e„_eoj 

Re7_^e„_e„jRe7+,<:„_eoj Re7+,e„_e„ jRe7+,<:„-eo j" 



j,mn 



Re7__ 



em-co,J 



Re7_ 



,«»i-«0 J 



Poo 



Re7_ 



Re7_ 



Re7_ 



(40) 
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which shows that, under the dynamical detailed balance condition, the intensities of the microscopic currents are 
uniquely determined by the single sequence Jj.mo- 

The following identities make the physical meaning of the currents Jj._mn{t) and Jfmni^) clear: 



Jm{t) : = 

d 



^ ^ I ^ ' Jj,mn{t) ^ ^ Jj,mn{t) 1 
' — 1,2 \n<m n>m / 



^^tr{\ern){ern\ps{t)] (41) 



is the difference between the quanta emitted from and absorbed by the level er. 



rn mj — l,2\n<m n>m / 

-ti-(^Hsps{t)) (42) 



d_ 

dt' 



expresses the fact that the variation of energy of the system is exactly balanced. 

On the other hand, the behavior of each microscopic current Jj^mn doesn't always follow a naive intuition. For 
example, even in the symmetric configuration of interaction (.gi(fc) = .92 (fc) = 9{k) and /ii — 112), there are cases when 
some micro currents flow backward (i.e. from the low to the high temperature reservoir), however it is impossible 
that all micro currents flow backward. A sufficient condition that the total energy current 



m \n<m n>m / 



4^^ (43) 



is positive when the reservoir 1 is at lower temperature than 2 is that 

^<1, Vm>n, (44) 

pnn 

i.e. that there is no inversely populated state. In addition, if all Ji,mn and J2,mn have opposite sign, the following 
strong relation [Gihbs domination hound) holds: 

Pnn 



However 



is not true when the stationary state of the system does not satisfy the detailed balance condition (See (p4D, (|25| ) and 
(|§)). In fact 

Jl,mn H~ J2,mn — ^ ^ Jj,mn 

i=i,2 

■/I'^mn + ■/2'^mn ^ 0. (47b) 

In other words, these stationary current can satisfy ( ^6| ) if and only if the stationary state of the system satisfies 
the detailed balance condition. When the stationary state can be described with detailed balance condition, the 
generalized temperature defined by (|2^) becomes constant which can be interpreted as the local temperature of the 
system in between two fields. Thus this condition gives a characterization of those non-equilibrium stationary states 
which are local equilibrium stationary states with current. We show an important example of such state in the 
following, however apart from few trivial cases, to satisfy the detailed balance condition strictly is impossible in this 
model as explained in the previous section. We consider the case where the detailed balance condition is satisfied 
approximately, i.e. the linear transport regime. 
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B. Linear approximation, local equilibrium and Onsager relation 



Here we show that the stationary current defined by (|3^) is consistent with weh-known non-cquihbrium physics 
in hnear regime. First we assume that the form factors in the interactions are the same for the two fields {gi{k) = 
92{k) '■= g{k)). This implies that the stationary solution is symmetric with respect to the indices 1 and 2. Now 
consider a small variation of these parameters 



Po = ^ , op ^ Pi- [32, and = ^ , 0^ = ^2-/^1 

and the first order expansion of the stationary solution in (5/3 and 5^. This gives 



(48) 



Pn 



Pi = A) + 2 ' A^i = Mo - ^ 



°ii I ^ 2 8/3i 2 a/32 2 dfii ^ 2 afj.2 



P2 = Po 



= Mo + 



Pi = Po, 111 = ^0 

P2 ^ Po, M2 = Mo (49) 



Using the symmetry (in 1,2) of pmm at 6/3 — Sp — 0: 
dpmm dp 

mm 



+higher order corrections 
dpmm dp 

n 



d(3i 



Pi = po, A*i = MO 
P2 = Po, M2 = MO 



9/32 



Pi = Po, Pi = A*o 

/?2 = /?0, M2 = MO 



9mi 



/3i = po, Ml = MO 
/32 = Po, M2 = MO 



dp2 



Pi ^ Po, Ml = MO ^^'^■^ 
P2 = /3o, M2 = MO 



Pi = /3o, Ml = Mo 
P2 = /3o, M2 = MO 



(corrections of order > 2). 



all the cross terms in ( [l9| ) cancel and we obtain 

Pmm o I 50 6fi Pmm 

Pi = Po + , pi = po — 

P2 = Po ^, M2 = Mo + "2^ 

Therefore as far as we consider J| up to the first order in 5p and dfi (linear transport regime) we can replace p„ 



(51) 



in the definition 



into 



Pmm — Pmm 



Pi = Po, Ml = Mo 
P2 = Po, M2 = MO 



Using 



Re7+,i,e 



■ Re7+,2,e„-£„ 
^mn 



- 5p4- 



= (em-Cn-A'o)/5o 



9/3o i9mo / e(""-^" -^0)^*0 - 1 

1 



7mn ( ^^^'^ 



d(3o dpoj e(^"-^"-^o)'^« - 1 
dk |go,e„-e„(fc)P^('^(fc) - (e™ - £«)) 



(52) 

higher order correction (53a) 
higher order correction (53b) 

(53c) 



fceo„ 



we get (we denote the approximate currents Jj^^) 



J2- 



2 

'fmri \ Pmm 
In 



J2. 



5(3 



d 



Sp- 



d 



g(em-e,i-Aio)0o 



1 



-^-(em - Cn - Mo) - OM 



d 



3(em-e„-A'o)/3o 



9(e, 
Mo) - (5m 



m - £«) e('--^"-W')/3o - 1 

d 1 

d{e,n ~ e„) e(^'"-^"-A'o)ft - 1 



^^2^1, mn 



(^m ^'0)^2 — ^l,mn 

In addition, if em ^ Po ^ Sfi, by the equilibrium approximation 

Pmm — " ^ , ^ — / ' 



(54) 
(55) 

(56a) 
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one can see 



j(E) _ J(E) _ _ j(E) , „ 

'-'l-^\,mn ~ ■-'l,mn ^ "^2.mn J- 



(56b) 



which hold the condition (|46^ 

From this i t is clear that, for the system S (say atom), the non-equilibrium effects appear as first order effects in 
the currents (56b), but only as second order terms in the state. This suggests a theoretical explanation of both the 
empirical success and the limitations of Kubo linear response theory. 

Let us show a relation ||l8| between the two currents Ji^2,mn and J^_^2 mn ~ '^i-^2,mn ~ tJ'oJi^2,mn, which is the 
analogue of the famous Onsager relation between the electric and heat currents in the conductivity problem. It is 
only an analogy because the carrier of our currents is a Boson particle and not Fermion (electron) . From (M) , 



Jl — >2,mn 








jQ 









6fio 

/3o 



(57a) 



where 



d 



mn \ t'mm 



Lr. 



' d{e,n - £«) e('= 

(e('=™ 



-Mo)/3o _ 1)2' 
n ran 



mm - Pnn) (< whcu ^ holds.) 



M„ 



or we obtain explicitly 



dJi. 



-(em 



d 



d6fi 



Mo)^r^ 



— L„ 



(57b) 
(57c) 

(58) 



which is the Onsager reciprocal relation. 

One can easily see that these currents produce positive entropy. Following ||l|, the entropy production with these 
currents is given as 



( 

V Po 



and as far as mm holds, since + TmnMmn = 0, 5S is positive for any (J/i, 5(3) except for 




6fi = 6/3 — 0, 6fi — (e. 



m 



Po 



(59) 



(60) 



which imply j[^2 r 



= 0. 



As is well known, Onsager reciprocal relation is understood as a consequence of microscopic symmetry of the 
dynamics, based on the following two assumptions p8| : (i) There exists an intermediate time scale between macro 
and micro dynamics, (ii) Average of spontaneous thermal fluctuation of the microscopic observable decaying is 
described by macroscopic transport theory. Notice that both the above assumptions were deduced in our model from 
the stochastic limit, (i) corresponds to the fact that the convergence to the stationary state of the system is described 
in the rescaled time scale. This time scale is exactly the time scale used in assumption (i). Moreover what the 
stochastic limit tells us is that the dynamics of the currents (or the transport coefficients) are given in terms of the 
time correlations of the original field in the initial state. This is nothing but the situation described by assumption 
(ii) . In the context of derivation of the Onsager relation between heat and electric currents by linear response theory, 
since there is no Hamiltonian which can describe the force generating a heat current whereas chemical potential can 
be treated always dynamically, (ii) has to be required as assumption[po| . In the present paper, both temperature and 
chemical potential are treated as parameters of the environment fields in the framework of the quantum mechanics 
for a open system. Moreover one should notice that the current is described directly in terms of the dynamics of 
the fields. It is also important to notice that the equilibrium state approximation (BO) is not necessary to derive 
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the Onsager relation (^8|). Usually, Onsager relation is derived assuming a symmetric property of the microscopic 
dynamics |p^ . However as is discussed in the next section, this symmetric property is equivalent to the requirement 
that the state is equilibrium (see below ([7l|)). Our results prove that the Onsager reciprocal relation (|8|) can be 
valid without any symmetry of the dynamics. Gabrielli, Jona-Lasinio and Landim illustrated such a possibility using 
a classical, solvable and phenomenological model plf . 

IV. MASTER EQUATION CHARACTERIZATION OF DYNAMICAL DETAILED BALANCE 

In the equilibrium case, it is well known that the detailed balance condition can be characterized by a generator of 
the master equation of the system interacting with the environment ||l^, . Given the dynamical semigroup which 
drives the state to an equilibrium state 



d 

where 



— Pt = C*pt, Pt Peq, (61) 



tv{XC*ipt)) = tT{ptC{X)). (62) 
The detailed balance condition or KMS condition for p^q is characterized by the following equations |l7) : 

tT{pegC+{A)B) :^ triply AC{B)) for all A, S (63a) 

C{X) - C+{X) = 2i[H,X] {H = H'^) for all X (63b) 

In this section, we prove a generalization of the above characterization to non-equilibrium stationary states in terms 
of the dynamical detailed balance condition defined in the previous section. 



We consider the forward and the backward Heisenberg evolution of a system operator X, i.e. (cf. |12| Chap I, 
section 1.1.29) 

jf\x) -.^UlxUt for t>0, j'l^\x) -.^U^tXUU for t<0 (64) 

where Ut is the time evolution operator in interaction picture. After stochastic limit and in the notations (17), (23), 
these lead to the master equations for observables 

Uj^iX))^ ^[A,(,f))] 



dt 



+r^+ (^-{E^El {j[^\x))} - E^{jf\x))El 
C{{3f\x))), for t>0 (65a) 



|(,f)(X))= z[A,(,f))] 



i^u- (\{ElE^,{j[''\x))} - El{j[''\x))E^ 



\{E^El {jf\x))} E^{j[^\x))E}\ 



-CB{{]f\X))), for i<0. (65b) 



where (•) denotes partial trace of the field degrees of freedom. Through (p2[), the dual master equation ( |22| ) (for 
density matrices) is written as 

j^ps{t)=C*ps{t), t>0. (66) 
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Similarly, we introduce a master equation associated to Cb as 



j/i\t) = -'^Bpf\t), t<0. (67) 



Both master equations have the same stationary state ps (see ( |23| ) and (^ 
As easily seen from (|65|), with A = given by (23) one has 

C{X) - Cb{X) ^2i[A,X]. (68a) 

By direct computation we obtain the deviation from the symmetry condition tr{psxC{y)) = tT{psCB{x)y) which 
characterizes equilibrium: 

tr {psXC{Y)) - tr {psCb{X)Y) = ^ XuY„,m (/9«(r_,e,-.„ + r+,,„_,j - p™™(r_,e„_,, + r+,,,_,„)) 

= ^'^XuYmrn(){ei — em)(JlJm + J2,lm) — ^(Cm — Q){Jl,ml + Jl.ml) (68b) 



where 
Choosing 
(p8q) becomes 



Xii = {<ii\X\ei) , Y^m ^ {tm\Y\tjn) , pii = {ei\ps\(ii). (68c) 
X ^\ea){ea\^.Pa, ^ = | Efc) (ef, | = : Pf,, (69) 
tr (psF,/:(Pb)) - tr {psCB{Pa)Pb) = 0{ea - eb){Ji^ab + J2^ab) - e{eb - ea){JiM + J2m)- (70) 



The left hand side describes the balance between two processes: transition from \ea) to jet,) and its converse in 
stationary state ps- Thus (|68| ) (or ([70|)) is a characterization of the dynamical detailed balance condition discussed 
the previous sectio n. R emember usual detailed balance condition is characterized by ( |63| ) which is the case when the 
right hand side of ( 68b ) is identically zero. 

Notice that ps is an equilibrium state when Ji.mn + J2,mn — 0. Let us remark again that as far as linear ap- 
proximation is concerned, Ji^mn = J2,mn = is not necessary to realize an equil ibrium state p^q (the equilibrium 
approximation (p9)) which follows the condition (p3) up to the first order (see Sec. [II B). In this case, 



tr {p,gXC{Y)) - tr {pegCB{X)Y) = tr {p,gXC{Y)) - tr {p,g£{X)Y) 

= (71) 

and it is exactly the symmetry of microscopic dynamics assumed in the original derivation of Onsager law[p^. 



V. LOCAL KMS CONDITION 



The KMS condition is known to be a characterization of equilibrium states equivalent to the detailed balance 
condition. In this section, we prove that a generalization of the KMS condition which characterizes the state described 
with the dynamical detailed balance condition. 

First, we introduce a generalization of the KMS condition which distinguishes between those general density matrices 
which commutes with a given discrete Hamiltonian and those which are function of the given Hamiltonian. This 
condition, which wc call local KMS condition in the sense of energy space, can describe states with mode-dependent 
temperatures 

Given a discrete spectrum Hamiltonian Hs'. 

Hs = J2eP, , P. = |e)(6| , Hs\e)^e\e) (72) 

For any complex valued Borel function / : M ^ C the map x i-^ g^^'tf{Hs)rj.^-itf{Hs) jg (defined by the spectral theorem 
and one has 

x{t) := f^^fiHs)^^-^tf{Hs) ^ J2 e^*(/(^)-/(^'»P,xP,, = ^ e'*'E{{x) (73a) 
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where 

%:={/(£) -/(£'); Ve,e'} , := ^ P^xP,,. 

e,e' : f(e)-f{e')=5 

For such Hamiltonian Hs the following theorem holds: 

Theorem 1. For a density matrix p and the corresponding state (( ■ )) the following are equivalent: 
(i) There exists a real valued Borel function /? : M ^ K. such that exp ~(3{Hs)Hs is trace class and 



(73b) 



^ Z 



-P{Hs)Hs 



(74a) 



(ii) There exists a real valued Borel function /3 ; R ^ R such that cxp ^P{Hs)Hs is trace class and p satisfies the 
following local KMS condition with respect to the Heisenherg dynamics x i— > e'*^^a;e~**^'^ ; 



Vx,y,i, {{xy{t + iP{Hs)))) = {{y{t)x)) 



(74b) 



where the meaning of y(t + iP(Hs)) is given by (73a) 
Proof. 

( |74i] ) => ( |74b| ). 



{{xy{t + tPiHs))) - tr (^pxe'^'-"^^"' y{t)e+^^"'^"-''^ = ^tr (^xe-^^"-'^"^y{t)'^ = tr iyit)xp) = {{y{t)x)) (75) 



( 74^ ) ( ^4aD . 

(74b|) means that for all x,y and for all t 



tr 



|^g-/3(ffs)Jfsy(^)g+/3(ffs)ffs^2,^ = tr(p2/(i):E) 



Therefore for all y and for all t 
or equivalently, putting t = and replacing y by ye~f^^^^^^^ 
hence, putting y = \ 
([78|) , ([79|) imply that, for all y 
and this implies that, for some scalar A 
Since tr(p) = 1, ( pT| ) implies that 



P~ Z 



(76) 

(77) 
(78) 
(79) 
(80) 
(81) 

(82) 



(Q.E.D) 

Notice that when j3{Hs) — (3 (constant), the state (74a) is the Gibbs state at temperature f3^^ and (74b) becomes 
the KMS condition. 

We shall prove that this local KMS condition (^^ is equivalent to the dynamical detailed balance condition (|6^). 
To avoid infinite-valued functions, we assume that all the pu are strictly positive and we represent the stationary 
solution ps of ( |66| ) and (|7|) in the form 



PS = le-'5-(^-)^^ psiei) = --logp„. 
Z ei 



(83) 



For such state the following theorem holds: 



16 



Theorem 2. The dynamical detailed balance condition ^6q) holds if and only if the local KMS condition ( [T^ j is 
satisfied. 

Proof (0) ^ (H). 

Appling the local KMS-condition ( [74|) to this state, we get 

{{AB)) ^ {{B{-^(3s{Hs))A)), (84) 

In addition in the notations (17), (23) and using relations 

Ai-il3siHs)) = A (85a) 



we obtain 

{{X[A,Y])) = {{XAY-A{-tf]siHs))XY)) 

= {{XAY-AXY)) 

= m,A]Y)) (86a) 

{{x{e1^^^^e,^^,^,y})) = {{xeI^_^^e,^_,y + xyeI^_^^e,^-,^)) 

= ((XEI^^^^E,,^^,,Y + El^^^J-iPs{Hs))E,,^^,J-il3siHs))XY)) 

= {{xe1^^^^e,^_,y + e1^^^^e,^.,^xy)) 

= {{{X,El^_,^E,^_,jY)) (86b) 

{{XeI^_^YE,^^,^)) = e''^(^")'"-^«(^'")^'"((Se„-e„X£;l^^_,^r)) (86c) 

{{XE,^_,YEI^^^J) = e'^'^''"'>''"-l^''^'-^'"{{El^_^^XE,,^_,^Y)). (86d) 

Now let us define Cq by the relation: 

{{C+{X)Y)) :^ {{XC{Y))) (87) 



for £ given by (65a). Notice that we are defining Cq not only in equilibrium state but also in the non-equilibrium 
stationary state which is described with the local KMS condition (Jt^), unlike ( |63a|) . Using relation (^, we find 

{uj — e,n — e„) 

c+{x) - ~i[A,x]~Y. (r-,- (^1{eIe^,x}-eIxeS^ +r+^. (^^{e^eIx} - e^xeI^^ 

+ J2 ((r+,^e'3s(^'")^"-^-(^")^" - T^,^)ElXE^ + (r^,,e'3-(^")^"-'3^(''")^" - r+,^)E^XEl) 
uieF 

= Cb{X) + nc.(^) (88a) 

ti^ix) = (r+,^e'3-(^")^"-''-(^")^" - r_,^)i?t xi?^ + (r_,^e'5-(^")^"-'5-(^-)^'" - r+, ji?^xii;t (88b) 

( p7| ) and (|8^) mean 

{{XL{Y))) = {{Lb{X)Y)) + J2 {{flUX)Y)) (89) 



uj£F 
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and 



/ p-0s(Hs)Hs , 



Cm -^Ti 

( sol) and (|9C) is exactly the dynamical detailed balance condition (|68|). 

Following (84)^(|90|) conversely, we see that the dynamical detailed balance condition (JG^) implies 

tr (psXCiY)) = tr (^pse^'^'^"''^"'' C{Y)e-^''-"'^^"'' , V X, F. (91) 
For off diagonal type operator 

there exists Y such that 

£(F) = y (92) 
and putting y — e**-'^sye~'*''^® (y is also off diagonal type) we get 

tr (p5^e"^^2/e-'*^-^) = tr [psei^siKsWs ^^iHs y^~^tHs ^-fisiiisWs ^ V 

or 

iT{psXy(i))=iT{psy(t^iiis{iis))X), VX (93) 

In addition, since p^t is diagonal, ( p3| ) is always satisfied with any diagonal type operator ?; = "Ylim ^ram\^m){^m\ also. 
Therefore, (^3|) is always satisfied with any operator X and y. 
(Q.E.D) 
Notice that since 

£+(l) = ^fl^(l)^0 

in the non-equilibrium case, cannot be a generator of any dynamical semigroup whereas Cb always exists as 
generator of dynamical semigroup. This is also one of the particular properties of the non-equilibrium state. In an 
equilibrium case, as we have seen Ps{x) become a constant /? which is the same inverse temperature of the environment 
fields, and the equality T^,i^/r^^i^ — e~'' *^'^™~'^"''holds, i.e. II^{X) = which implies Cb — Cq. 

VI. INTERACTION WITH NON-EQUILIBRIUM FIELD 

In the previous sections, we considered the non-equilibrium stationary states of a system driven by two environments 
at two different temperatures and we discussed several characterizations of such states. In this section, applying these 
characterizations to the state of the environment, we consider a system interacting with an environment in local 
equilibrium. (On the local KMS condition for the field degrees of freedom, see the next section.) One will see not only 
that the stationary state of the system driven by such non-equilibrium environment can be characterized as for the 
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previous model, but also that interesting non-linear effects due to the interaction with non-equilibrium environment 
exist whose physical meaning is different from the previous model. 

We consider a system interacting with a single boson field whose state is described by a generalized temperature 
function. Technically, the analysis of the model can be done in the same way as the previous one. Instead of (jlj) but 
similarly, the Hamiltonian 

H = Hq + XHj, (A is a coupling constant.) (94a) 



Ho = Hs + Hb, Hs^}^ 

I 



ei\ei){ei\, Hb = J u;ik)aiak [at, al] ^ Sik ~ k'), (94b) 



H, = dk 



(^g{k)Dal+g*{k)D^ak) . (94c) 



On the other hand, we assume that the initial state of the field is a mean zero gauge invariant Gaussian state with 
correlations: 

{aia,,)^Nik)6ik-k% Ar(fc) = (95) 

where (3{u){k)) is some positive function. This is a natural generalization of the Gibbs factor to which it reduces when 
P{uj) is constant: 

fi{Lo) = p. (96) 
Exactly in the same way as in the previous argument, one can derive the white noise Hamiltonian equation 

where 

bt:u.= [ dk g:{k)bt:Jk), 5t..(/c)= limie-*('^^ W-")*/^'afe. (98) 
The state of the limit white nose will be of the same type with correlations 

(6L(fc)&t':-'(fc')> = S^^-27rS{t-t')d{k-k')6iLo{k)~to)N{k) (99a) 
{btUk)bl,,^,{k')) = 5^^,27rS{t~t')5ik-k')S{Lo{k)~Lu){N{k) + l). (99b) 

Finally we obtain the master equation ( ^2| ) but with different parameters 

A = t J2 {lMl-u.)ElE^ - lm{^+^)E^El) (100a) 
Fzp,^ = 2Re 7zF,^ > 0, (Fz^,^ = for u;<0). (100b) 

where 



= / dk\g^{k)\ 



2-i{N{k) + l) 



Lo — U!{k) — iO 

= / dk\g^{k)\ 



,2 -^N{k) 



bj — uj{k) — iO 

TT [ dk \gjk)\^——l— 5(uj(k) - w) - iP.P [ dk . (lOOd) 
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As in the previous model, the ofF-diagonal elements vanish when Gmn 7^ 0, (V m,n) which is defined in (|23c|). In 
order to see if the stationary state can violate the detailed balance condition or not, let us check condition (^61). With 
direct computation we find 

^ g+/3(e^-.,)(^™-^0 for em>ei (101a) 

1 Im 



El!^ ^ for ,^<,, (101b) 

1 Im 

Let us remark this fraction does not depend on the structure function g{k) unlike the previous model, however it can 
violate condition ( p6| ) due to the generalized temperature function /3(w), i.e. 

Tm/ / Tmfc ^ kl /ino^ 
Im ^ km ^ Ik 

except for the constant temperature case (^6|). 

Let us show a typical example of non-equilibrium effects due to the generalized temperature function. To realize the 
stationary state with non-detailed balance condition at least, two Bohr frequencies (three level system) are necessary. 
With a generic 3-level system, whose energy levels are given by ei < 62 < 63 and 63 — 62 7^ 62 — 61: the concrete form 
of the matrix A in ( p^ ) is written as 

r+,£2-£i + r+^gg^gj — r_^e2-ei ~r_.e3_.(:j 

and one can directly see that its eigenvalues are 



A = 0, ^ > 0, (104a) 

b = r+_E2-ei + r+,e3_ei -I- r+^e3_e2 -I- r_^e2_£j -I- r_^g3_ej -f r_^e3_£2 (104b) 

C = r+,e2-ei-'^+,C3 — £2 -'^+,C2 — ci ,C3 — ei + -'^+,<S2 — ci ,C3 — £2 ^~ -'^+,C3 — ei ,C2— ei "I" -^+,£3 — ei -'^+,£3 — ^2 

-t-F-i- £3_£jr_ £3_£2 -|- r_ — eir_^e3_£j -|-r_ ^2— eir_^e3_e2 "I" T-l-jea— £2-^ — ,£3 — ei (104c) 



and the stationary state 



where 
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Pll 
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— ei 
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r 
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-ei 


^r_,£3 


p 

— ei-^ +,C3 


-£2 


P33 


r+,e3- 


P 

£2-^ +,£2 


— ti 


+ r+,£3- 


F 

Cl-^ -,f2 


-ci 


f r+,£3 


p 

-£2-^ +,£3- 


-ei 


P22 


r-,£3- 


■eir+,e2 




-f r_^£3_ 






^-r+,£3 


P 


-£2 



(105a) 

X (105b) 

y (105c) 

Z. (105d) 

When ( ^6[ ) is not satisfied, the above solution does not satisfy the detailed balance condition. Notice that in this case 
the detailed balance condition is equivalent to 

8 /3(62 - 6l)(62 - 61) - /?(63 - 6l)(63 - 61) + - 62)(63 - 62) = 0. (106) 

Let us remark that the physics of this model can be different from the previous model. For example, taking (6i|Z?|62) = 
(so as r±^£2_£i = 0) for simplicity, the above quotients become 



P22 _ r+,£3_£ir_ 



Pll r-,£3-eir 



■■<^3-t2 _ g/3(£3-£2)(£3-£2)g-/3(£3-ei)(e3-ei) X (107a) 



,£3 — ei-*- +,e3 — «2 
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P33 ^ T^+,ea^ci ^ ^^fj(e3->ii)ie3->^l) ^ P33 ^ ^ +.£3-62 ^ ^-I3(e3-e2){e3-e2) ^ ^ (107b) 

pii r_c3_gj P22 r_(:3_e2 



and X is larger than 1 when 



/?(e3 - £2) > ^^^/3(e3 - ei) (108) 
£3 - £2 



Thus, for such temperature function /3(x) the stationary state satisfies P22 > Pii which means that 2 is a so-called 
inversely populated state. 

Here, we focus on the current passing through the stationary state and discuss the non-linear effects. For simplicity, 
we discuss the case of a three level system. In this case (ei < £2 < £3), with direct computation we obtain 

Jmn ■ = / dk{U}nuUt) 

/3(e2-ei)(e2-ei)-/3(e3-«l)(e3-ei)-l-/3(e3-e2)(e3-e2) _ 1 

_ (_'l\m+n+l ^ T (^^^Qa\ 

^ > (e/5(E2-ei)(e2-ei) _ ;^)(g/3(e3-£2)(e3-£2) _ _ e-/3(e3-ei)(e3-ei))^ ^'''''"^> 

Jin- ^ I dkL0{k){U}nkUt) ^{t^-en)J^n (109b) 
I = \{e,\D\^2)? I dk\g{k)\H{u{k)~{e2-e,)) 

X \{e2\D\e3)\' [ dk\g{k)\^6{u;{k)-{e3-e2)) 

x\{e3\D\e,)\^ f dk \g{k)\H{oj{k) ~ {es - ei)) (109c) 

Notice Jjf ^ and Jg^''' have same (and Jgf ^ has opposite) sign. In addition 

4f^ = -i4f^ + 42^) (110) 

and the sign of each currents depends on 

S /3(e2 - ei)(e2 - ei) - /3(e3 - ei)(e3 - d) + /3(e3 - e2)(e3 - £2). (HI) 

In the case 5 = 0, all currents vanish. Especially, when the function is a constant (3 (i.e. the initial state of the field 
is an equilibrium state with temperature this is easily understood with the fact that the state of the system 

converges to the equilibrium state at the same temperature without any stationary currents. Notice that even within 
the linear approximation up to order S, there is no local stationary state (with currents) which satisfies the detailed 
balance condition, unlike the previous model. In this model, the existence of currents always implies the deviation 
from the equilibrium. 

Now let us see some interesting properties of the currents (109). In the case S > 0, we obtain from (109) 

4f\-^if >0, and j}f <0. (112) 



As clearly understood from the definition of the currents, the relation ( |112| ) is describing the process that a field 
quantum with energy €3 — ei is converted into two quanta with energy €2 — ei and €3 — 62. On the contrary when 
S<0, 

4i\42^<0, and Ji^f >0 (113) 

and this can be interpreted as a process from two quanta to one quantum. There are interesting analogies of these 
processes with parametric downconversion and second harmonic generation in non-linear quantum optics |^2j. They 
are considered as opposite process of another. In our model, the direction of the process depends on the generalized 
temperature function P{uj) which is a parameter of the initial state of the field. This phenomenon can be understood 
as the fact that through interaction with a non-equilibrium field the system can have such a function, which is an 
example of dissipative structure in the Prigogine sense jl] . 
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VII. DISCUSSION 



In conclusion, let us further comment on a few related topics. 
1) On the irreversibility and unitarity of time evolution. 

due 
ing the 



As we discussed in Sec. [II B, we can see irreversibility in this model through the entropy production 
to the stationary currents, which should be considered as processes involving the total system includi 
environment. On the other hand, the time evolution operator of the total system Ut is unitary in the sense that 



uJUt = UtUj - 1, t > 0, (114) 

which is easily checked by putting X = 1 in (pT|). These statements might seem to be contradiction. However, 
one should notice that the appearance of irreversibility has nothing to do with the unitarity of Ut . When the 



temperatures of both environments are the same, it is known that the unitarity condition (114) is required 
to realize a physical fluctuation-dissipation relation or a correct equilibrium stationary state |12 . Moreover, 
when we speak of macroscopic phenomena like entropy production, we need a good procedure to extract the 
proper degrees of freedom to discuss them. Since there exist same macroscopic states which are distinguishable 
microscopically from each other, not all microscopic degrees can be employed to discuss macroscopic properties. 
Indeed, the entropy production (|5|) is discussed in terms of what we call slow degrees of freedom, and the 
stochastic limit can be considered as the procedure of extracting the proper degrees of freedom. In other words, 
we extract information from the total dynamics as slow degrees which can describe the macroscopic phenomena. 

2.) Local KMS condition for field. 

A possible formulation of the local KMS condition for the field is the following. 

Definition 1. A state (•) on the polynomial algebra au, a\., is said to satisfy the local KMS condition with 
temperature function f3 : — > M if for every m, n G N, Ei, . . . , e„, r]i, . . . , rjj^ € {0, 1}, and with the convention 
x'^ = x'^ , — X for any operator x, the following identities hold in the sense of distributions. 



K (0) • • • alZ {Q)a'C + '(3hJall-_\ {t + */3^„. J . . . {t + i^, )) 



lHt)...al\{t)al(Q)...alim- (115) 



Lemma 1. Define the local inverse temperature function by 

n{i 

m{k) J ujk ' 



/3(fc):= llog^^)-, (116) 



where 

{aual,) =: m{k)5{k - k') = = (gn(A:) + \)8{k - k') (117) 



(q ~ —1 for Bosons and q = +1 for Fermions). Then the local KMS condition is satisfied by the 2-point 
functions: 

{ak{0)al{t + tP{k')))^{al{t)ak) (118) 

{ai{0)ak'{t + z/3(fc'))> = K'(t)4(0)) (119) 
Proof . In the above notations, one has 

(afc(0)4,(i + z/3(fc'))) = e^(*+^^('='))-'=' (a^at,) = g-z^C^')-.' e^*-.' (a.al) 



e 



'"^''^n{k)S{k~k') = e'*^-'m{k')S{k' ~k)^e'*^-'{alak) = (4,(i)«fe) 



and this proves (118). In a similar way one verifies that (|119|) holds 
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Proposition 1. // the state (•) is mean zero gauge invariant and Boson Gaussian then condition (Hi) is 
satisfied. 



Proof. By Gaus sian ity both sides of ( 115 ) are reduced to weighted sums of pair correlation functions. Since 
in both sides of (115) we can distinguish the (/i, e)~terms from the (fc,77)-terms and since the pair correlations 
preserve the order, there will be 3 types of pair correlations: (i) those of type {h, k), (ii) those of type {h, h), 
(iii) those of type {k,k). 

In case (i), due to gauge invariance, the only none zero combinations are of the form (aa^) or (a^a) so we can 



apply dTTa) and (119) 



In case (ii) the terms are already in the correct order. 

In case (iii), again by gauge invariance, the only possibilities are 



= {al{t)ah>{t)) (120) 

and similarly for the other term. 

Since in the Boson case the weight of each pair partitio n is equal to 1, after the replacements (|11^), (119), (12C) 



the pair-partition expansion of the left hand side of (115) becomes the pair-partition expansion of the right 
hand side. 

The validity of the local KMS condition for more general Gaussian states as well as for quantum Markov states 
is now under investigation. 

3.) The generalized temperature function and its thermodynamics. 

On the description of the generalized temperature function P{H), R. S. Ingarden, A. Kossakowski, M. Ohya, 
T. Nakagomi have discussed similar idea in the context of information theory They introduce a system 

described by the density operator 

and discussed possible generalization of thermodynamics for structured complex systems (e.g. biological system) 
including bifurcations, catastrophes and self organization. As mentioned in their book 123], their phenomenolog- 
ical idea is in the line of thought of synergetics by Haken |24| . In the present paper, we explained the microscopic 
origin of such states and their physical meaning through the dynamical detailed balance condition. Through the 
local KMS condition a general classification of such non-equilibrium states became possible. We believe that 
our approach gives a good insight to generalization of thermodynamics in this direction. 
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